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■r::j- ' Abstract 
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c^ \ We use the Thomas-Fermi method to examine the thermodynamics of particles 

. , obeying Haldane exclusion statistics. Specifically, we study Calogero- Sutherland par- 

Q \ tides placed in a given external potential in one dimension. For the case of a simple 

O ■ 

harmonic potential (constant density of states), we obtain the exact one-particle spa- 



/\ • tial density and a closed form for the equation of state at finite temperature, which are 

both new results. We then solve the problem of particles in a x'^'^ potential (linear 
density of states) and show that Bose-Einstein condensation does not occur for any 
statistics other than bosons. 



PACS numbers: 5.30.-d, 71.45.Jp 



A recent definition of generalized exclusion statistics by Haldane has aroused con- 
siderable interest [1]. The definition is based on the rate at which the number of avail- 
able states in a system of fixed size decreases as more and more particles are added to 
it. We introduce a statistics parameter ot and assign the values a = and 1 to bosons 
and fermions respectively, because the addition of one particle reduces the number of 
available states by a. Some examples of particles obeying such generalized statistics 
are spinous in an antiferromagnetic spin chain with inverse square exchange [1,2] and 
two-dimensional anyons residing in the lowest Landau level in a strong magnetic field 
[3-5]. 

More recently, the statistical distribution function for such systems has been ob- 
tained [5,6]. If [I is the chemical potential , e is the energy of a single-particle state and 
T is the temperature, then n(e, /i, T) satisfies 

(n-^ - aT in-^ + 1 - af'"" = exp[(e - fi)/T]. (1) 

(We set Planck's and Boltzmann's constants equal to one). Let g{e) be the density of 
states so that the total number of particles and the energy are given by 

N = de g{e) n{e) 

\ (2) 

and E = de g{e) n(e) e. 

Using Eq. (1), the thermodynamic potential n(/i,T) satisfying (9n/9/i)T = — A^ can 
be shown to be given by 

n = -T I de g{e) log[l + exp(- i/T) ], (3) 

where e(e,/x, T) is defined by the relation 

e + T (1 - a) log[l + exp(- i/T) ] = e - /x. (4) 

In terms of e, n is given by 

1 

^ = ,~ /r^X • (5) 

exp(e/r) + a ^ ^ 
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If we keep T fixed and let /x ^ — oo, we find that e ^ e — /x and n(e) -^ exp[(/i — 

e)/T] independent of a. In this hmit, therefore, A^, ^ and E all go to zero as exp(/i/T). 

We now use the identity 

dn dn ^ ,^, 

& + s^I = «• («' 

which follows from Eq. (1) if T is held fixed. Now suppose that g{e) ^ e" . Inserting 
(6) in (2), we obtain 

(f ), = (" + 1' ^- (^' 

Hence E + {a + l)Q is independent of /x. Since £" + (a + 1)0 is zero at /x = — oo, it 
must be zero for all /x and T. This relation will be used below. 

The simplest example of particles obeying the generalized exclusion statistics is 
provided by the Calogero-Sutherland (CS) model [7,8]. This one- dimensional model 
has the Hamiltonian 

ff = J_ V P? + ^^^^^^^ V , ' „ . (8) 

2m ^^-^ m ^^-^ (xi - XjY 

Here a > with the understanding that the wave function vanishes as jx^ — Xj |" 
whenever the particles i and j approach each other. Note that a = and 1 yield free 
particles (bosons and fermions repectively) while the maximally attractive potential 
occurs for a = 1/2 (semions). Using the thermodynamic Bethe ansatz [9] and other 
means, it has been argued that the CS system can be thought of as an ideal gas in the 
sense that the particles only have statistical interactions amongst each other [6,10,11]. 
While the CS model has been studied in great detail over the years, we would like 
to examine it here from the view point of exclusion statistics. The CS model has so far 
been exactly solved only in an external simple harmonic potential V{x) = mto'^x'^ /2 . 
However, once we realise that the particles only undergo statistical interactions amongst 
each other, we can develop other techniques to study the thermodynamic of the CS 
system placed in any external potential . One such technique is the Thomas-Fermi 



(TF) method which has been used to study fermionic systems for a long time [12]. 
In addition to being exact in the large-A^ hmit, the TF method does not require a 
knowledge of the energy spectrum and wave function of the system. Here we modify the 
method by using the distribution in (1) rather than the Fermi-Dirac distribution while 
continuing to use the usual phase space description of quantum systems. We show below 
that this procedure leads to the well-known expressions for the one-particle spatial 
density and the ground state energy for the many-body system in a simple harmonic 
potential at zero temperature [7,8]. Encouraged by this, we analyse the situation at 
finite temperaure and obtain the one-particle density which is a new result. We also 
find closed forms for the thermodynamic quantities /x and Q, in contrast to previous 
expressions in the literature which are given as high-temperature series. Finally, to 
demonstrate the power of the TF method, we study the CS model placed in a x"^'^ 
potential which has g{e) ^ e. Once again, we derive the thermodynamic quantities 
for all a and T. Since the g{e) is the same as for free particles in four dimensions, it 
is natural ask whether Bose-Einstein (BE) condensation occurs at low temperature for 
a finite range of a. It is shown that condensation occurs only for bosons (a = 0) and 
not for any positive a. 

The TF density of states is obtained from the integral 

dx e{e - V{x)) ^^' 



n 



^2{t - V{x))lm 



where 6'(y) = 1 if ?/ > and if y < 0. The TF expression for the one-particle density 
is 

p(x,r) = I ^n(e,T), (10) 

where e = p^/2m + V{x) and / dx p{x,T) = N. 
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At T = 0, Eq. (1) implies that there is a Fermi energy /x(0) such that 

n{e) = l/a if e < /x(0), and 

(11) 
= if e > /i(0). 

Thus a < 1 (or > 1) corresponds to less (or more) exclusion than fermions. We now 
study the harmonic potential problem. From (9), g{e) = l/a; . Therefore the Fermi 
and ground state energies are given by 

/i(0) = aojN 

1 (12) 

and ^(0) = - iV/i(0). 

Eqs. (10, 11) yield the well-known semicircle law [8,13] 

p(x,0) = — (2m/i(0) - m^iJx^fl'^ (13) 

for x^ < Xq = 1[i{^) I muJ^ . The thermodynamic limit for this system is obtained 
by taking N ^ oo holding /i(0) and therefore u>N fixed. (This is analogous to keeping 
the density in a box fixed with l/a; playing the role of the volume). 
At finite temperature, we find from Eqs. (2, 6) that 

where n(0) is the solution of (1) with e = 0. This has the solution 

fi{N,T) = aujN + T log[l - exp(-a;A^/T)]. (15) 

(Here we used the boundary condition that for fixed T, A^ ^ as /x ^ — oo). Now 

^'"^ Ni^] (16) 



implies that [14] 



dN J J. \dNyj. 



1 rp2 pujN/T 

n = -aojN^ + — dy ^. (17) 
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Since E = —Q, the entropy S (given hy 0, = E — TS — fiN ) is independent of a. Eq. 
(17) has a high-temperature expansion in u>N/T of the form 

-n = NT + ^LoN^ {a - ^) + ... . (18) 

This resembles the virial expansion for the equation of state in a box and it has the 
remarkable property that only the second term depends on a [11]. However, this 
expansion only has a finite radius of convergence (27r). Eq. (17) has an alternative 
expansion in powers of exp(— iVw/T ) which is valid for all T. For instance, we may 
read off the low-temperature result 

-n = \aujN'' + ^— (19) 

z do; 

plus exponentially small terms. 

We may now obtain the finite temperature density p{x,T) from Eqs. (1) and (10). 

This is somewhat difficult to compute explicitly except for special values of a. For 

instance, 

2 

n(e) = — -^^=^^ for a = 1/2 and 
V 1 + 472 

1 1 r^r~ ^^^ 

= ^ ~ ^\ -. for a = 2, 

2 2V4+7 

where 7 = exp[(e — lj)/T] . In Fig. 1, we show p{x,T) at T = ujN/2 for a = 1/2, 1 

and 2. For convenience, we have rescaled x -^ x/xo and p -^ p{x,T)/p{0,0) so that 

p{x,0) = Vl-x^ 9{l-x^) 

(21) 



and 

'0 



/ dx p{x,T) = 7r/4 

for all a. Although T is fixed in Fig. 1, its ratio to p(0) = atoN is greater for smaller 
a. This explains why the curve for a = 1/2 is spread out to the largest extent. We 
have also exhibited the T = (or a -^ 00) semicircle for comparison. 

We now apply the TF idea to a CS system placed in a different potential which has 
not been studied before. An interesting and exactly solvable TF case is the potential 

V{x) = i (9mr]^x^] . (22) 



The TF density of states is given by g{e) = e/rj . Hence E = —2Q. At T = 0, we 

find that 

/x(0) = (2aT]N 



1/2 



and ^(0) = - Nfi{0). 



(23) 



The thermodynamic hmit is again defined by taking N ^ oo keeping r]N fixed. To 
study finite temperature, we define an auxihary variable 



oo 



i/ = / de n{e). (24) 

Jo 



Eq. (6) imphes that 



/'dN\ 
\dfxjj, 

and TT— = n(0). 



(25) 



As in Eqs. (15, 17), we deduce that 

lJ,{iy,T) = au + T log[l — exp(— i//T)] 

1 r/'^ y (26) 

and riN{v,T) = - au^ + T^ dy ^ _ ^ . 

Using (7), we can compute the energy 

Eiu,T) = 2 / dy N{y,T) ^ ii{y,T). (27) 

Jo oy 

Thus /x, A^ and E are all known in terms of f and T. 

Since a particle moving in a two-dimensional harmonic oscillator potential or mov- 
ing freely in four dimensions also has a linear density of states, we may ask whether our 
system exhibits BE condensation at low temperature for a finite range of a. If a = 0, 
we see that (26) has no solution once T falls below 



T, = ^/G^tTV/tt. (28) 

For T < Tc , u and /x stay at oo and respectively. Therefore the second equation in 
(26) can only be satisfied if there is a macroscopic occupation A^(0) of the e = state, 
that is, if r]N = r]N{0) + tt'^T^/G . 
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However, for any a > 0, Eq. (26) has a solution right down to T = and there is 
no condensation at finite temperature. In Fig. 2, we present /i vs T with tjN = 1 and 
a = 0, 0.02, 0.1 and 0.25. It will be seen that n{T) is nonanalytic only for a = at 
T = v^/tt = 0.7797... . 

To summarize, we have developed a TF method to study the thermodynamics of 
particles obeying the generalized statistics of Eq. (1) in any given external potential . 
We can numerically compute the one-particle density for any a and T. If g{e) varies 
as a (nonnegative) integer power of e, we can solve the thermodynamics exactly. For a 
linear density of states, there is no BE condensation if a > 0. This also implies that a 
1/x^ interaction between bosons, no matter how weak, destroys BE condensation. 

One of us (D. S.) thanks the Physics and Astronomy Department of McMaster 
University for its hospitality during the course of this work. This research was sup- 
ported by a grant from NSERC (Canada). 
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Figure Captions 

1. The one-particle density p{x,T) vs x in a harmonic oscillator potential with T = 
ioN/2 and a = 1/2, 1 and 2. The T = (or a ^ oo ) semicircle is shown as a 
solid line. 

2. The chemical potential /x vs T for a linear density of states with r]N = 1 and 
a = 0, 0.02, 0.1 and 0.25. The solid line (a = 0) is nonanalytic at T = \/6/7r. 
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